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Метод вычисления функций 
в неавтономном режиме 


Предложен таблично-алг оритмический мегод вычисления функций с использованием смещенных избыточных 
систем счисления, который позволяет совмещать процессы поразрядного ввода аргумента и поразрядного 
формирования результата. Это дает возможность использования метода для неавтономного выполнения 
зависимых по данным операций в режиме совмещения. Показана зависимость времени вычислений от 
параметров систем счисления. 


Введение 


При реализации последовательно-параллельных алгоритмов в параллельных вы- 
числительных системах можно выделить последовательность (цепочку) операций (функ- 
ций), которые нельзя распараллелить ввиду их зависимости по данным. Следующую 
операцию в цепочке можно начать выполнять только после получения результата пре- 
дыдущей операции. Аналогичная ситуация возникает при реализации итерационных 
вычислительных процессов, когда на очередном шаге вычислений используется резуль- 
тат, по крайней мере, одного предыдущего шага. 

Для ускорения вычислений в данном случае может быть использован метод неав- 
тономного вычисления цепочки зависимых операций в режиме частичного совмещения. 
Зависимые операции, принадлежащие критическому пуги на графе алгоритма, выпол- 
няются в цепочке операционных устройств, между которыми данные передаются по- 
разрядно [1]. В каждом устройстве реализован метод выполнения операции, позволяющий 
на каждом шаге совмещать во времени ввод операндов, их обработку и выдачу одного 
разряда результата в следующее устройство. При таком режиме вычислений следующая 
операция начинает выполняться сразу после получения первого разряда предыдущей 
операции. За счет этого обеспечивается частичное совмещение зависимых операций, 
что создает предпосылки для ускорения вычислений [2]. 

Современные методы проектирования вычислительных систем на основе ПЛИС 
дают возможность на одной микросхеме создавать композицию устройств переработки 
информации, позволяющую выполнять одновременно операции над многими 
операндами. Эффективность использования ресурсов ПЛИС во многом зависит от 
способов обмена данными между компонентами системы. Последовательная передача 
данных между компонентами сокращает необходимое число внутренних и внешних 
связей ПЛИС, что при прочих равных условиях позволяет улучшить характеристики 
использования ресурсов ПЛИС. 

Наличие в составе ПЛИС блоков памяти дает возможность для ускорения вычисле- 
ния функциональных зависимостей использовать таблицы функций. Известны методы 
неавтономного воспроизведения функций с использованием таблиц для симметричных 
систем счисления с цифрами {-п,...,-1,0,1,...ли} [3]. Однако наличие цифр с отрицатель- 
ными знаками приводит к расширению разрядной сетки для представления промежу- 
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точных результатов и к увеличению объема аппаратуры, причем как для устройства 
воспроизведения функций, так и для других устройств, работающих совместно с данным 
устройством для неавтономного выполнения цепочки операций, зависимых по данным [4]. 
Устранение указанных недостатков является актуальной задачей, решение которой позво- 
лит повысить эффективность неавтономных вычислений. 

Ниже предлагается таблично-алгоритмический метод, который обеспечивает 
воспроизведение функций в смещенных системах счисления с различными основаниями. 


Обоснование метода 


Пусть функция У’ = Г(Х) задана в первом квадранте и является непрерывной. 


Функция может иметь экстремумы в области определения. Для определенности будем 
считать, что аргумент и функция представлены в избыточной смещенной системе счис- 
ления в виде: 


х=УхА- е Ук, (1) 
1 


где ху, Е {0,т} - цифры; А- основание системы счисления. Соотношение А < т 
определяет избыточность системы. 
Учитывая, что цифры результата формируются после обработки входных дан- 
ных, то есть с запаздыванием, будем рассматривать функцию У =А ”/(Х). Здесь р- 
задержка в числе шагов формирования разрядов результата относительно ввода соот- 
ветствующих разрядов аргумента. 
Естественно, что для получения и разрядов после запятой функции У не- 
обходимо сформировать т = п+ р разрядов функции У. 
Значение функции, зависящее только от 7 старших разрядов аргумента, обозначим 
ар 
ВЕКУ, (2) 
где Х, - код аргумента, представленный 1 старшими разрядами. 
Если на каждом шаге вычислений значение У, выбирать из условия 
-р 7 
КТ а, (3) 
то на п-м шаге погрешность вычисления будет меньше А”. 
Введя обозначение 
_ р-р — 
В =" У(Х)-УК, (4) 
условие (3) можно записать в виде 
0< № <1. (5) 
С учетом (2) и (4) можно заключить, что перед началом вычислений (1=0) 


условие (5) выполняется. Предположим, что оно выполняется на 1-м шаге. Опреде- 
лим, как обеспечить выполнение на следующем шаге условия 


0< А <1. (6) 


1-1 


Для (1+1) -го шага с учетом (1) и (4) можно записать 
В ЕЕ (7) 
Значение функции /(Х,+х„.К"") представим в виде Г(Х,)+А 


1 
= К’ ‚то есть 


м 
через значение функции на предыдущем шаге и приращение, отнесенное к разряду с 
весом К "'. Тогда (7) преобразуется к виду 


ВАНА Аа (8) 
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Введя обозначение 


Ни = ЕВ, + "Ад, (9) 
условие (6) с учетом (8) запишем как 
ОН, <1+Ущ. (10) 


Левое соотношение в (10) выполняется, поскольку функция рассматривается в 
первом квадранте. Для выполнения правого неравенства необходимо обеспечить вы- 
полнение соотношения 
т-к+1 

А 


тах 


ть : (11) 
которое получено с учетом (5) и (9). 

Выражение (11) дает возможность определить значение задержки р, если из- 
вестно максимальное значение приращения функции А„„ в области ее определения, 


или определить диапазон изменения цифр избыточной системы, то есть максимальное 
значение т при заданном значении р. 

Таким образом, подбором значений р и т можно обеспечить выполнение усло- 
вия (3) на каждом шаге вычисления функции, что гарантирует получение результата с 
заданной погрешностью. 

Несложный анализ формулы (11) показывает, что система счисления обязатель- 
но должна быть избыточной (т > ^ ), а задержка формирования разрядов результата не 
может быть меньше единицы. Заметим, что увеличение диапазона изменения значений 
цифр нежелательно, поскольку это приводит к увеличению разрядности шины для пере- 
дачи цифр между устройствами. 


Формирование таблицы приращений 


Когда формирование приращений функции на каждом шаге нельзя вычислить 
аналитическим путем или это приводит к большим затратам времени, целесообразно 
использовать таблицу приращений, хранящуюся в памяти. 

Таблица формируется следующим образом. 

1. Строится ярусно-параллельный граф перехода значений аргумента, состоящий из 
п-+1 яруса (включая 0-й). Каждая вершина графа 1-го яруса соединяется дугами с т+1 


вершиной (1+1)-го яруса в соответствии со значениями цифр аргумента х, е {0 т}, 
которые приписываются указанным дугам. Каждая вершина имеет обозначение С’ ‚ где 
1 (0, п) номер яруса, а / е {0,К' — !— номер вершины на 1-м ярусе. 

2. Каждой вершине п -го яруса присваивается вес С’, = "У (5%. 


3. Для ярусов с номерами 1 < и веса вершин определяются по формулам 


. м 
ПС ии |" — 0,1...., ти) ыы (п-1) (и-2) 0 
т , , 20. 


> > 1 у м 
4. Для каждой дуги графа, соединяющей вершину С’ с вершиной С,’ , опре- 


С, = еп 


деляется приращение 
7 р м 7 
А(С',, С )=(, -КО,. 
Приращениям 1-го яруса придается вес А". 
Покажем, что сумма всех приращений, принадлежащих дугам каждого пути от началь- 


ной вершины до каждой конечной, равна значению функции для соответствующего 
значения аргумента в избыточном представлении. 
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Пусть путь проходит через вершины С”, С", С2,...,С”",С”. Тогда сумму 5 прира- 

щений с учетом их весов можно записать в виде 
1 Оу -1 2 1 —2 п п-1 —п 
5 =(С.-АС, Ко +((, - КО +. + (С, -КС” К”. 

Раскрывая скобки и учитывая, что С) =0, получим 5 = (”К" =У(Х,), то есть 
сумма приращений дает действительное значение функции для заданного значения аргу- 
мента. 

Адрес приращения в таблице однозначно определяется дугой графа, а дуга, в 
свою очередь, определяется совокупностью верхнего и нижнего индексов вершины 
а, ‚ для которой дуга является входящей. 


Алгоритм и пример вычисления функции 


Для аппаратной реализации целесообразно использовать основание системы, 
равное 2* (2 = 1,2,3...). В этом случае могут использоваться обычные двоичные аппа- 
ратные компоненты (сумматоры, регистры и т.д.), а умножение на величину основа- 
ния заменяется сдвигом. 

Полученные табличные значения приращений могут превышать значение макси- 
мально допустимой цифры т избыточной системы счисления. Поэтому алгоритм вы- 
числений предусматривает преобразование суммы приращений в последовательность 
цифр у, Е {0,т}. 

Процесс вычисления функции на каждом {1-м шаге определяется формулами (8) и (9). 

В соответствии с формулой (9) находится Н,. Целая часть этой величины яв- 
ляется очередной цифрой результата, то есть у, =епё Н,, а сдвинутая влево на один 
разряд дробная часть используется на следующем шаге в качестве ^,. Начальным яв- 
ляется значение А, =0. 


Пусть необходимо вычислить непрерывную функцию У’ = /(Х) ‚ заданную в пер- 
вом квадранте, с использованием смещенной избыточной системы счисления с основа- 
нием К =4 и цифрами ху, Е {0,4}. Для представления аргумента ограничимся и =3 
разрядами. 

Полный граф в данном случае имеет 5 вершин на 1-м ярусе, 25 вершин на 2-м 
ярусе и 125 вершин на 3-м ярусе. Для примера рассмотрим фрагмент графа, который 
определяется частью таблицы значений (табл. 1). 


Таблица 1 — Фрагмент таблицы задания функции 


Значение аргумента Значение функции 
к=2 К=4 К=2 К=4 
0,110000 0,234 0,111111 0,333 
0,101111 0,233 0,111000 0,320 
0,101110 0,232 0,110000 0,300 
0,101101 0,231 0,111011 0,323 
0,101100 0,230 0,111100 0,330 


На рис. 1 показан фрагмент графа, полученный в соответствии с приведенной выше 
процедурой формирования таблицы приращений. Нижняя дуга каждого разветвления 
соответствует цифре аргумента х, = 0, а верхняя — х, =4 (чтобы не загромождать рису- 
нок, цифры аргумента приписаны только двум дугам первого яруса). Приращения, соот- 
ветствующие дугам, показаны жирным шрифтом. 
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0-й ярус 


6‘ =0 


2-й ярус 


62 =20 || 


6% =20 | 


62 =12 | 


62 =10 [ 


Рисунок 1 — Фрагмент графа 


3-й ярус 

Со, =333 
С =320 
С =300 
0: =328 
С. =330 
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Для полученного фрагмента графа максимальное приращение в десятичном пред- 
ставлении равно Л. =15. По формуле (11) определяем минимальную задержку фор- 


мирования цифр результата р =2. Для получения значения функции У’ = Г(Х) необ- 
ходимо выполнить р+п = 5 шагов. 


Для значения аргумента Х = 234 процесс вычислений иллюстрируется табл. 2. 


Таблица 2 — Состояния переменных 


1 х, А, Н(=2) | Н(=4) | у Микрооперация 
0 =: _ 00,0000 0,00 _ _ 
+00.0001 +0.01 суммирование 
1 2 1 00,0001 0,01 0 
00,0100 0,10 сдвиг 
+00,1000 +0.20 суммирование 
2 3 20 00,1100 0,30 0 
11,0000 3,00 сдвиг 
+00,1111 +0,33 суммирование 
3 4 33 11,1111 3,33 3 
11,1100 3,30 сдвиг 
+00.0000 +0.00 суммирование 
4 - — 11,1100 3,30 3 
11,0000 3,00 сдвиг 
+00.0000 +0.00 суммирование 
5 в _ 11,0000 3,00 Э 
00,0000 0,00 сдвиг 


Целая часть Н,, показанная жирным шрифтом, является очередной цифрой 


результата. После выполнения 5 шагов получено 5 дробных разрядов функции 


У=К°Г(Х) = 0,00333. Начиная с третьего шага формируются разряды функции 


У = Г(Х) = 0,333 . Полученный результат соответствует указанному в табл. 1 значению 
функции для заданного аргумента. 
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Заключение 


Предложенный таблично-алгоритмический метод воспроизведения функций по- 
зволяет совмещать процессы поразрядного ввода аргумента и поразрядного формирова- 
ния результата в смещенной избыточной системе счисления. Это дает возможность 
использовать метод для неавтономного выполнения зависимых операций в режиме сов- 
мещения, когда другие операции в цепочке выполняются также в смещенной системе 
счисления. 

По сравнению с методами, использующими симметричные системы счисления, в 
данном случае упрощается аппаратная реализация и уменьшается время вычислений, 
поскольку нет необходимости работать с отрицательными приращениями. На практике 
широко применяются вычислительные методы, например, численного интегрирования, 
цифровой обработки сигналов, вычисления полиномов, которые в ряде случаев могут 
быть реализованы в смещенных системах счисления. 

Благодаря поразрядной передаче информации уменьшается число связей между 
компонентами системы по сравнению с использованием методов параллельной арифме- 
тики. Это особенно важно при реализации систем на ПЛИС, поскольку повышает надеж- 
ность систем и более экономично использует ресурсы интегральных схем. 
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ТА. Дичка, В.В. Жабйна 

Метод обчислення функщй у неавтономному режим! 

Запропоновано таблично-алгоритмачний метод обчислення функщй з використанням змицених надлишкових 
систем числення, що дозволяе сумпцення процесв порозрядного введення аргументу 1 порозрядного 
формування результату. Це дае можливсть використання методу для неавтономного виконання залежних за 
даними операшй у режим! сумпцення. Показано залежесть часу обчислень в1д параметрав систем числення. 


Т.А. Русйка, Г.Т. йабта 
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